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1 Introduction 

One of the main reasons to introduce the curvaton scenario [1-4] was to obtain a reahstic 
mechanism of generation of non-gaussian adiabatic perturbations of metric [1, 5]. Since that 
time, many interesting curvaton models were proposed. However, it would be nice to have 
a curvaton model which would be as simple as the basic chaotic inflation scenario with the 
potential m?'(jP' /2 [6]. It would be good also to find a natural implementation of this scenario 
in the context of supergravity. This is the main goal of our work. 

In this paper we will describe a broad family of models of chaotic inflation in super- 
gravity, which provide a natural realization of the curvaton theory. We will calculate the 
non-gaussianity parameter /nl for the simplest versions of these models and show that this 
parameter takes different values in different parts of the universe, in agreement with the cur- 
vaton web scenario of Ref. [7] . If inflation is sufficiently long, the average value of /nl in this 
scenario does not depend on the initial value of the curvaton field. We will also show that 
under certain conditions the parameter /nl can take values in the observationally interesting 
range from 0(10) to O(IOO). 
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2 Curvaton scenario and chaotic inflation in supergravity 



For many years, it seemed very difficult to realize chaotic inflation in supergravity. This 
problem was solved in [8]. The authors proposed a very simple model describing two fields, 
S and with the superpotential 

W = mS^, (2.1) 

and Kahler potential 

/C = 55-^(^-^)2. ^2.2) 

Note that the Kahler potential does not depend on the phase of the field S and on the real 
part of the field <1>. Therefore it will be convenient for us to represent the fields S and <I> as 
S = (Te*^/\/2 and ^ = {(p + ix)/V^- The field (p plays the role of the infiaton field, with the 
quadratic potential, as in the simplest version of the chaotic infiation scenario [6]: 



2 

V{())) = 3H^ = ^</)2 , (2.3) 



where H is the Hubble constant during inflation. Near the inflationary trajectory with 5 = 0, 
the mass squared of the imaginary part of the field $ is = GH'^+m?. Thus during inflation 
> 6H^, and therefore the imaginary part of the fleld <^ is stabilized at Im<l> = 0. No 
perturbations of this field are generated. 

Both components of the field 5 may remain light during inflation, and therefore in- 
flationary perturbations of these fields can be generated [9]. Since the potential does not 
depend on the field 6, we will ignore fiuctuations of this field in our study of the curvaton 
perturbations. The potential of the fields o" at x = is 



.2 

V{cP,a) = —e 



m ^2/2 



(2.4) 



For a <^ 1 one has 



m<x) = ^ + ^ + ^. (2.5) 

The effective mass squared of the field cr at a ^ 1 is given by 

2 2,3 2 i2 2 2,9 2 2 /r, n\ 

m„ = Vo-cr = m + -m (p a = m + -H a , (2.6) 

where V^o- means second partial derivative of V with respect to a. One can easily see that 
= for (pa <^ 1. During inflation <C H^, and therefore inflationary perturbations 
of the fleld a can be generated. At (pa > 1, the effective mass squared of the fleld a is 
dominated by the term ^m'^(p'^a'^ = ^H'^a^ > m^. For cj ^ 1 one still has ^ H^, so 
the perturbations of the fleld a are generated in this regime as well. However, at cr > 1 the 
potential becomes exponentially steep, and ^> H^. Therefore inflationary fluctuations of 
this field are generated only for u < 1. This is a very important advantage of the curvaton 
scenario in supergravity: The steepness of the curvaton potential at a > 1 protects us from 
extremely large perturbations of the curvaton field which otherwise could be produced during 
eternal inflation in this scenario [7, 10]. 
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If one does not take into account the curvaton fluctuations in this scenario and study only 
the usual inflaton fluctuations [11], then the COBE normalization requires m ~ 6 x 10~^, in 
the system of units Mp = 1 [6, 12-14]. Thus, the mass of the inflaton field must be somewhat 
smaller than 6 x 10^^ if we want to add the curvaton fluctuations to the inflaton fluctuations. 

Recently the supergravity model described above was substantially generalized in [15, 
16]. The generalized scenario describes a theory with a superpotential 

W = Sf{<^), (2.7) 

where f{^) is an arbitrary real holomorphic function. The Kahler potential in this class of 
models may take several different functional forms, e.g. 

JC = SS-^{^-^)^-^iSSf. (2.8) 
In this theory, the inflaton potential is given by 

V{^) = fWV2). (2.9) 

and the mass of the field a is 

ml = aH' + {f{HV2)f. (2.10) 

In this class of models, one can implement chaotic inflation in supergravity, with an arbitrary 
shape of the inflaton potential V{(j)). In all of these models one has = f{^/V2)/3. The 
term (/'(0/\/2))^ is equal to 3H^e, where e <C 1 is the slow roll parameter. For a> 1 one has 
> H^. In this case no curvaton perturbations are produced, so all standard predictions 
of the single- field inflaton scenario remain intact. 

On the other hand, in models with a <^ 1 one has <^ during inflation, which 
means that quantum fluctuations of the field a are generated during inflation [15, 16]. 

Thus we have a broad class of models of chaotic inflation where the curvaton scenario 
can be realized. One can further generalize this scenario by adding terms ~ to the 
superpotential, and by using other versions of the Kahler potential, as long as the Kahler 
potential has certain properties described in [15, 16]. The requirements which are necessary 
for the existence of the light curvaton flelds in this class of models can be formulated in an 
invariant way in terms of the curvature of the Kahler geometry. In particular, the parameter 
a is related to the curvature of the Kahler manifold [16]. The fleld a itself has an interesting 
interpretation from the point of view of supergravity: it is the scalar component a of the 
goldstino multiplet. Because of the generality and simplicity of this scenario and because of 
its supergravity origin, one may call it the supercurvaton scenario. 

In this paper we will concentrate on the simplest model (2.1), (2.2), but with an addi- 
tional term —■^{SS)'^ in the Kahler potential, as in Eq. (2.8). In this model the curvaton 
mass squared along the inflationary trajectory with cj = is given by 

ml = m'^ + aH'^ , (2.11) 

and in a more general case < cj <C 1 the effective mass squared of the fleld a is 

ml = m^ + aH^ + ^H^a^ = + ^m^^ + ^m^cj^^a^ . (2.12) 
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3 Curvaton perturbations and non-gaussianity 



During inflation, the curvaton perturbations are produced. An average amplitude of pertur- 
bations produced during each Hubble time is given hy 6a = Then these fluctuations 
are stretched, overlap with each other, and eventually produce a classical curvaton field a 
which look relatively homogeneous in the observable part of the universe, but may take dif- 
ferent values in other parts of the universe [7] . The amplitude of the perturbations of density 
of the curvaton field with a quadratic potential is given by 5pa/pa ~ 16a ja. However, the 
total energy density of matter at the moment when the curvaton field decays may be greater 
than the energy of the classical field a. This may happen, for example, if the decay of the 
inflaton field during reheating produces many curvaton particles [7]. Therefore the relative 
perturbation of density will be given by 

' > (3-1) 

P 0" 

where r = Paj p at the time of the curvaton decay. According to [5], these perturbations will 
match the COBE normalization of the spectrum for 

r— ~ 7 X 10"^ . (3.2) 
a 

These perturbations are non-gaussian, with the amplitude of local non-gaussianity given by 
[5] 

/nl = |: . (3.3) 

Our goal will be to find a typical value of a in some of the simplest supergravity models 
described above, calculate Jo", find the value of r required to satisfy Eq. (3.2), and finally 
determine /nl- The most complicated part of this program is finding a typical value of a. 



4 Stochastic approach 

We will begin our study with investigation of the behavior of the distribution of the fluctua- 
tions of the curvaton field a with a simple quadratic potential m'^a'^ /2. This approach will 
allow us to describe the case when = m? + aH^ , but not the more general situation when 
depends on a as in (2.12), which will be discussed separately. 

During inflation, the long-wavelength distribution of this field generated at the early 
stages of inflation behaves as a nearly homogeneous classical field, which satisfies the equation 

ma + V„ = Q. (4.1) 

or, equivalently. 



However, each time interval H 
average amplitude squared 



da^ = . (4.2) 

dt 3H ^ ' 

new fluctuations of the scalar field are generated, with an 
i^-') = S • (4-3) 



-4- 



The wavelength of these fluctuations is rapidly stretched by inflation. This effect increases 
the average value of the squared of the classical field o" in a process similar to the Brownian 
motion. As a result, the square of the field a at any given point with an account taken of 
inflationary fluctuations changes, in average, with the speed which differs from the predictions 
of the classical equation of motion by ^ 



dt 3H 47r2 

Using 3H<j) = —V^, one can rewrite this equation as 



+ (4.4) 



It's solution with the initial condition a{(f)i) = for = m? + aH^ is given by 

/ 



(4.5) 



where (hi is the initial value of the inflaton field. 



2 



F'(0) 



(4.6) 



If inflation continued for much longer time than 60 e-foldings, as we will assume in 
this paper, the main contribution to a is given by perturbations produced at the very early 
stages of inflation. Such fluctuations look almost absolutely homogeneous on the scale of the 
observable part of our universe, so our calculations give us a typical value of the classical 
field a inside the observable part of our universe. However, in different parts of the universe, 
the field a may be significantly smaller or greater than its "typical" value calculated above. 
As a result, the amplitude of the curvaton perturbations is not a constant, but it varies in 
space [7]. The same is true for the degree of non-gaussianity /nl, see Section 6. 

Therefore, to be precise, one should distinguish between the average amplitude of the 
field a calculated above, when the averaging it taken all over the universe, and the local 
value of the field a in each horizon-size part of the universe. We will make this distinction in 
Section 6, where we will make a slight change of notation and call the value of the curvaton 
field averaged over the whole universe a, reserving the letter a for the average value of the 
curvaton field in the horizon-size part of the universe. However, in the main part of our 
paper we will not distinguish between a and a. This means, in particular, that when we will 
calculate /nl(c), our results will in fact describe the value of this parameter for a = a, i.e. 
the value of /nl for an average value of o", all over the universe. In Section 6 we will show, 
however, that the value of /nl for an average value of a can be significantly different from 
the average value of /nl; the order of averaging in certain cases can be very important. One 
should take this effect into account when making predictions of the non-gaussianity in each 
particular curvaton scenario. 

We should note also that in general the curvaton field may not be equal to zero at the 
beginning of chaotic inflation, so one may also consider a possibility that initially {a'^){(pi) 
was very large. In this respect, the supergravity model which we are going to study provides 
an important simplification: The curvaton field initially cannot be much larger than 0(1) 
because of the exponential steepness of the potential at a > 1. Also, the effective mass 
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term |m^(/)^(T^ in the supergravity potential (2.4), rapidly reduces the initial value of the 
field (T, thus making quantum fluctuations generated during inflation more important than 
the initial value of the classical field a. We will study these issues in the rest of the paper, 
starting from the simple toy model with = and ending up with the model with 
ml = m^ + ^m'^(fP' + |m^(/>^(T^. As we will see, in all these cases the final result does not 
depend on the initial distribution of the curvaton field if infiation lasts long enough. 



5 Nongaussianity in various regimes 

5.1 A toy model with = m? 

In this subsection we will study the distribution of the curvaton field with the mass = m^ 
during infiation driven by the massive inflaton field with potential V = ^m^cj)^. In this case 
Eq. (4.6) implies that the classical scalar field a which is nearly homogeneous on the scale 
of the horizon has a typical amplitude 

Meanwhile the amplitude of fiuctuations of a generated at that time is 

^ H m(j) ,^ 

During the subsequent evolution of the universe, a and Sa both decrease in the same way, 
and therefore at the end of inflation the curvaton perturbations have fiat spectrum with the 
amplitude 

- = I- (^•') 

a (pi 

As we mentioned above the amplitude of the perturbations must be normalized as 

^ ^ ^ 7 X 10-5 . (5.4) 

a (pi 

and hence 

5 3.5 X 10^ 

/nl = ~ (5.5) 

This means that the degree of nongaussianity depends on the initial value of the inflaton 
fleld. Unless this fleld is very large, /nl may be extremely large. 

However, in supergravity models which we study in this paper the approach developed 
above is valid only if inflation was short enough, that is, (pi <ti m~^/'^ and for large values of 
(pi one cannot ignore the supergravity correction to the mass in Eq. (2.6). 

5.2 m2 = m2 + fj/V^ 

In the previous section we made a simplifying assumption that the curvaton mass does not 
depend on a, which allowed us to use Eq. (4.6). However, as one can see from (2.4), in 
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supergravity model (2.1), (2.2) the curvaton mass does depend on a in a rather comphcated 



way. The leading correction to the curvaton mass squared m? is given by 'jm'^cfP'a'^ = %H^a'^ 



and it becomes dominant for (f)a > 1. 

To find out how it will change the final result one has to solve Eq. (4.5) for 

which takes in this case the following form 

y = - + ^ - , (5.7) 

X 4 

where x = (j)'^ , y = a'^ , and b = 

The general solution of this equation can be expressed in terms of Airy functions, 

^(^)- (2i,)2/K Ai(^)-cBi(z) 

y(x)- [2b) ^Ai'(z)-cBi'(z) ' ^^-^^ 

where z = 2~'^^^b^^^x and c is a small constant which should be chosen in such a way that 
y{xi) = 0. 

Suppose first that the initial value of the field (p is much higher than m^^^'^, i.e. z ^ 1. 
One can check that in this case one should take c ^ 1 to have y{xi) = 0. Infiation ends at </> ~ 
1, when z <^ 1. In this limit, all functions are 0(1). Therefore the functions Bi(z) drop out 
from the final expression because of the small coefficient c, Ai(z) ^ Ai(0) = S'^/^ r-^(2/3), 
Ai'{z) ^ Ai'-i/3r-i(l/3). As a result, 

(25)V3^r(l/3) 

3i/3r(2/3) • ^'-'^ 

Expressing everything in terms of the original fields (f) and a, we find 

rj,-. ^^^^^ r(l/3) 2/3, 

a(6) « — ^ 0.15 m^''^ 6 . 

24/3V37r2/3 r(2/3) ^ 



(5.10) 



This yields 

— ~0.4m^/^. (5.11) 
a 

The COBE normalization requires rm^/^ j x 10^^. Therefore, for m ~ 10^^ we have 
r ~ 0.04 and /nl ~ 30. 

If, on the other hand, the initial value of field 4> is much smaller than m~^^^, then the 
final result looses it's universality and become sensitive to In this case, one can either use 
the analytical solution above, with different initial conditions, or simply use the results of 
the previous section (one can see that in this case (pa <^ 1, and hence the results of Section 
5.1 are valid). 

Note that in our calculation of /nl we used equation (3.3), which was obtained in [5] 
under the assumption that the curvaton potential is purely quadratic. Meanwhile in our 
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case the curvaton potential contains the quadratic term as weU as the quartic term 

2)2 4 

{'g°^ , see (5.6). This could lead to some corrections to equation (3.3) [17]. Fortunately, 
one can show that during the last 60 e-foldings of inflation in our model the quartic term is 
vanishingly small as compared to the quadratic term. That is why one can use the simple 
equation (3.3) for the calculation of /nl- 



5.3 = + cxH^, a > 



Now we will consider the case when the mass of the curvaton field is given by 



ml = aH"^ =m^ + 1^ , (5.12) 

where we have ignored the correction ^rn^ cjP' a'^ , which will be taken into account in Section 
5.4. 

To study this case we consider separately the evolution of perturbations at acj)^ /Q > 1 
and acfP' < 1 assuming that during last 60 e- folds of inflation the condition acjP' /Q < 1 is 
satisfied, which means that a < 1/40. Thus during last 60 e- folds, ~ m?, and hence one 
can use the results of Section 4. 

Substituting m? = aH^ = in (4.6), we obtain 

/ 



^ 1 /• 1/2(0) 



v 



2a f V 



(5.13) 



For the case of the power-law potential V the integral in (5.13) can be calculated exactly. In 



particular, for V = \w?clP' and a (j),"^ ^ m cfp' one obtains 



/ 2 6 



ibiT^a \ a J 

Note that that this result does not depend on the initial value of the inflaton field. At the 
end of the first stage of inflation when a4>i/Q = 1, both terms in the brackets are equal to 
each other and the averaged value of a at that time is about 

, , , v3m 

^(•/-i) ^ ^ , (5.15) 
while the amphtude of the perturbations of field a is 

777 

5a{ct>i) ^ . (5.16) 



The CMB normalization of the amplitude of the perturbations thus requires 

r— ~r./f ~7x 10-5 . (5.17) 
a V 3 

Note that in this case the amplitude of the curvaton perturbations does not depend on the 
inflaton mass m. 
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Taking a = 10 ^ we find that r ~ 10 ^ and hence /nl ~ 10^. Meanwhile for a = 10 ^ 
we should have r ~ 10~^, which gives /nl ~ 10^. 

One may wonder what is the origin of such an incredible sensitivity of the results to the 
choice of the parameter a. The answer is that this parameter makes the mass of the curvaton 
field much greater than the mass of the inflaton field at the early stages of inflation. As a 
result, the distribution of the field a shrinks fast while the field (p rolls down. 

In the calculations above we have ignored the supergravity correction to the curvaton 
mass squared: 'gm?(j)'^a'^ = ^H'^a^. As we will show in the next section this correction can 
be ignored only if a ^ 10~^m^/^ and hence the results of this section are applicable only in 
this case. 

5.4 ml = m^ + cxH^ + ^H^a^, ct > 

Now we will study the curvaton perturbations in the theory with the general potential 
which corresponds to the curvaton mass (2.12). For o"^ ^ 1 equation (4.5) becomes 

f^y + -^ + y^.,^, (5.19) 

ax X 6 4 



where x = cfP' , y = cr^, b = In this case, unlike to Eq. (5.7), there is no exact analytical 



967r2 

solution. Nevertheless one can investigate the solutions of this equation using phase diagram 
method. If infiation lasts long enough, all solutions, independently of the initial conditions, 
converge at a certain attractor trajectory in the phase space {y,x), or, equivalently in the 
space {a, (j)), see Fig. 1. For large (p, this attractor trajectory is given in the leading order by 
the solution of the algebraic equation: 




which is 

The existence of the attractor solution implies that if inflation is long enough, the final 
results do not depend on the choice of initial conditions for the curvaton field. We have also 
found above that in the limit a — )■ one should get the asymptotic solution (5.10), whereas 
for large a the asymptotic solution is given by (5.14), (5.15). One may wonder how large 
should a become for the switch between these asymptotic regimes? 

To answer this question, let us use the variables: 

x = z 6-^/3 , y = u h^/^ , a = 7 h^'^ , (5.22) 

in terms of which equation (5.19) becomes 

dz z ?, A 
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Figure 1. Behavior of the average value of the curvaton field cr as a function of the inflaton field 0, 
for various initial conditions. As we see, all trajectories which start at the early stages of inflation 
(large field cj)) converge to the same attractor solution. We follow it until the field 4> becomes 0(1) 
and inflation ends. At large 4>, this solution is very close to the square root of the function (5.21), 
which is shown by the blue dashed line. 

After rewriting equation (5.19) in this form it becomes clear that the behavior of the solutions 
is controlled by a single parameter 7 = ab~^f^ ~ 10 a m^^/^. One can easily understand 
that the two asymptotic regimes discussed above corresponds to 7 ^ 1 and 7^1. One can 
confirm this conclusion by direct numerical calculations. 

This means that the results obtained in Section 5.2 are valid for a <^ 10~^m^^^. Mean- 
while in the opposite limit a ^ lO^^m^/'^ one should use the results of Section 5.3. To give 
a particular example, let us take m ~ 10^^ In this case one can use the results of Section 
5.2 for a <^ 10~^, whereas for a S> 10~^ one should use the results of Section 5.3. 

5.5 ml = m^ + aH^ + ^H^a"^, cx < 

Finally, we will study the case a < 0. At first glance in this model the mass squared of 
the curvaton field at large and cr = is negative, and therefore one expects a tachyonic 
instability. However, similar to the model considered in the previous section, one can show 
that for \a\ ^ 10~^m^/^ the effect related to the negative mass squared contribution aH^ is 
subdominant and can be ignored. In this case the results obtained in previous Section 5.2 
are applicable. 

For lO^^m?^^ ^ \a\ ^ 1, the tachyonic instability leads to spontaneous symmetry 
breaking controlled by the supergravity correction ^H^a^ = '^m?(jP'a'^ to the curvaton mass 
squared. Indeed, one can show that the minimum of the supergravity potential for the 
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curvaton field, in the regime with |a|,cT <C 1, can be found from the following equation: 

9 2|a| 4 

<y =^ + 13 ■ (5.24) 



Therefore at large (j) and a < the potential has a minimum at 

= ^ . (5.25) 

This means that during inflation the field a falls towards this minimum, and its distribution 
become centered not at o" = but at o" = As for the height of the potential along the 

trajectory with a = y'^^^' ^"^^ small |a| it remains approximately given by 12. 

When the field <^ becomes smaller than 6/|a|, the minimum of the potential shifts 
towards o" = 0, and the curvaton mass squared becomes equal to m} . However, this does 
not mean that the distribution of the field a instantly follows the position of the minimum. 
Since the mass of the field a at that time is much smaller than iif, the field a will move 
towards <t = very slowly, decreasing at the same rate as the amplitude of perturbations 
bo. As before, we are assuming that |a| < 1/40, and therefore the curvaton mass squared is 
given by rr? during the last 60 e-folds of inflation. This leads to the following result for the 
perturbations: 

bo 

— ~ — ^ , , . (5.26) 
o 2^/27r|a| 

For ra ~ lO"'^ and a ~ — 10 a proper amplitude of perturbations corresponds to r ~ 1/3 
and, hence, /nl = 0(3). However, one can easily increase /nl by increasing m and/or 
decreasing For example, taking m ~ 10^^ and a ~ — 10~^ gives /nl ~ 30. 



6 Non-gaussianity and the curvaton web 

In the previous sections we have evaluated the average value of the curvaton field at the 
last stages of inflation, and calculated the parameter /nl describing local non-gaussianity. 
However, we should remember that when we were talking about the classical homogeneous 
curvaton field cr, we had in mind the long-wavelength perturbations which look homogeneous 
on the scales corresponding to the present observable part of the universe. In reality this 
classical field in our model is a random variable with the expectation value o = \J {bo"^) 
obtained by summing up the contributions of all long wavelength fluctuations (larger that 
the present horizon) generated on inflation. All calculations above were performed taking o 
to be equal o. However, because cr is a random gaussian variable it take different values in 
different parts of the universe of the size of our horizon [7] . 

To evaluate the observational implications of this fact, let us try to understand how the 
amplitude of perturbations of metric and the local value of /nl depend on the local value of 
o. For simplicity, we will assume that the standard inflaton perturbations are very small, so 
that we can ignore them in our investigation. This can be achieved by considering a model 
with m ^ 6 X 10~^. We will also assume that the curvaton field density at the moment of the 
curvaton decay is much smaller than the total density, i.e. r <C 1. In this case, the change 
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of a does not affect the total density p, but it does affect Sp{a), which is proportional to a. 
This means that the amplitude of perturbations of metric produced by fluctuations of the 
curvaton field will be proportional to |: 

p p a ' 

Meanwhile the local value of /nl is inversely proportional to r = For small a, the value 
of p(cr) is proportional to o"^. Therefore 

-2 

/nl(^) = /nl(ct) • ^ . (6.2) 

The probability that the curvaton field will take some value much greater than a, is 
exponentially small. However, the probability that a is substantially smaller than a can be 
rather large. To estimate this probability we will make a simplifying assumption. Namely, 
we assume that all values of the field S = with l^l < a are equally probable, but 

the probability vanishes for |5| > a. In the maximal value of the curvaton field is (J\/2, the 
probability to find the field a in the interval da from to is given by and the 
average value of the curvaton field is a, as it should be. 

Now let us evaluate the average value of the amplitude of density perturbations, averaged 
over all possible values of a: 

/W \ ^ ^faada ^ 2^2 5p(a) ^^^^^ 



\ p / p J a a"^ 3 p 



Thus, the average amplitude of the curvaton perturbations almost exactly coincides with the 
amplitude of perturbations in the universe with an average curvaton field a. 

However, the situation with (/nl) is quite different. Since /nl is proportional to cr"^, 
it expectation value over the whole universe acquires a divergent contribution from the parts 
of the universe with small a. Our calculations are valid only for fluctuations produced well 
before the last 60 e- folds of inflation, with a combined amplitude cr above 0{H). Introducing 
the cut-off at fj ~ ~ 2TT6a, we find 

(M ^ (.) / (^)' ^ = /NL (.) ... (^) . (6.4) 

H 



How significant is the effect discussed above? To give a particular numerical example, 
let us consider the case a = 0. In this case, according to Eq. (5.11), one has ^ ~ 0.4 
We found that for m ~ 10^'' one has /NL(f^) ~ 30. In this case Eq. (6.4) implies that 

(/NL)~5/NL(fT)~150. 

Thus we deal with a significant effect of statistical amplification of non-gaussianity: 
Even though the fraction of the volume of the universe with /nl(<7) ^ ftihi^^) is relatively 
small, the values of /nl in those parts of the universe can be huge, so the expectation value of 
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/nl can be much greater than the value of this parameter fNL{(^) calculated in the previous 
sections of the paper. 

This effect becomes even stronger in the models where the curvaton field is real (instead 
of being a radial part of a complex field) . In such models 

(/nl) ^ /nl (^) [(-Y^^ /nl (ct) (6.5) 

H 

In the particular example discussed above, ^ ~ 0.4 m^/^ and m ~ 10^'', this would 
lead to an enormously large amplification effect: (/nl) ~ 10^ /NL(f'') ~ 3000. 

Thus we see that in the curvaton scenario some fraction of the universe can be in a 
state with the curvaton field a significantly smaller than its average value a. In such parts 
of the universe, the locally observed level of non-gaussianity will strongly exceed its value 
/NL(f^) calculated in the previous sections. This effect is so significant that the average value 
of the parameter /nl can be much greater than the value /nl in the part of the universes 
with an average value of the field a. In other words, operations of averaging in this case are 
not commutative. 

For a complete investigation of this effect one should also take into account the standard 
inflationary perturbations of metric. The curvaton perturbations are important only in the 
cases where the standard inflaton perturbations are suppressed. That is why we assumed 
that m ^ 6 X 10~^. But the standard inflaton perturbations may dominate in the rare parts 
of the universe where cr <C ex. In such cases one should perform a more detailed investigation 
of nongaussianity of perturbations produced by all sources. 

This means that one should be very careful when formulating predictions for the non- 
gaussianity parameter /nl in the curvaton scenario, because the distribution of possible 
values of /nl in the curvaton web can be extremely broad. Moreover, the existence of the 

anti-correlation between the amplitude of the perturbations of metric {^^'^^ and the non- 
gaussianity parameter /nl for r ^ 1 (see equations (6.1) and (6.2)) suggests that anthropic 
considerations may play a very important role in evaluation of the probability to live and 
make observations in parts of the curvaton web with different values of the non-gaussianity 
parameter /nl [1, 7, 18-21]. We hope to return to the discussion of this issue in a separate 
publication. 

The difference between (/nl) and /nl(^) clearly demonstrates that /nl is not a perfect 
tool for the description of non-gaussianity. As shown in [7], the distribution of the regions 
of small (large) perturbations of metric and spikes of non-gaussianity has an interesting 
structure, which we called "the curvaton web." This structure has a non-perturbative origin. 

Indeed, the non-gaussianity parameter /nl(o") takes its largest values in the regions of 
the universe where the classical curvaton field a is small, see (6.2). In the theories where the 
curvaton field is a real, single component field, the regions of small a correspond to domain 
walls separating large domains with cr > from large domains with o" < [7]. 

In the theory studied in the present paper, the curvaton field a corresponds to the 
radial component of a complex field S. In this case, the regions of small a form strings. 
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reminiscent of the cosmic strings which appear due to spontaneous symmetry breaking. In 
our case, however, unlike in the usual cosmic string case, the curvaton strings appear in 
the places corresponding to the minimum of energy of the curvaton field. If one considers 
more complicated models, where the curvaton has 0(3) symmetry, instead of the domain 
walls and cosmic strings one will have localized objects reminiscent of global monopoles. In 
other words, the distribution of the peaks of non-gaussianity in the curvaton scenario has 
topological origin, which cannot be fully described by the standard tools of perturbation 
theory, such as /nl and (7nl- 

7 Discussion 

In this paper we discussed the curvaton scenario, which naturally emerges in the simplest 
supergravity realization of the chaotic inflation scenario [8, 15, 16]. Investigation of this sce- 
nario consists of several parts. The main step is to find an average value of the curvaton field 
a after a long stage of inflation. One needs this to calculate the amplitude of perturbations 
of density of the curvaton field. We performed this investigation by analyzing the growth of 
the curvaton perturbations during inflation. 

To conclude this investigation, one should find the ratio r of the energy of the curvaton 
field to the energy density of all other particles and fields at the time of the curvaton decay. 
This is a complicated and model-dependent problem, which requires study of reheating after 
inflation, the decay rate of the curvaton field, and the composition of matter at the time of 
the curvaton decay. In this paper, we simply treated r as a free phenomenological parameter, 
but one should remember that all of the issues mentioned above should be addressed in a 
more detailed investigation. 

We analyzed the model with the simplest quadratic inflaton potential and with the 
curvaton mass given by aH^ + m^. Our investigation demonstrates that if infiation is long 
enough, then the average value of the curvaton contribution to the amplitude of metric 
perturbations, as well as the averaged value of the non-gaussianity parameter /nl, do not 
depend on initial conditions for the curvaton field. The final results depend on the inflaton 
mass m, and on the parameter a, which is related to the curvature of the Kahler manifold 
[16]. However, the locally observable parameter /nl and the amplitude of the curvaton 
perturbations may take different values in different parts of the universe and in certain cases 
they may significantly deviate from their averaged values [7]. Moreover, the average value 
of the parameter /nl can be much greater than the value /nl in the part of the universes 
with an average value of the field a. For a certain choice of parameters, the value of the 
non-gaussianity parameter /nl can be in the observationally interesting range from 0(10) to 
O(IOO). 

The curvaton perturbations in our simple model have flat spectrum. This is a conse- 
quence of degeneracy of the masses of the inflaton and curvaton field at the end of infiation. 
One can change the spectral index by switching to a theory with a different infiaton potential. 
This can be easily realized in the new class of chaotic infiation models in supergravity, or by 
splitting the spectrum of fluctuations of the curvaton fleld into two branches with different 
masses [15, 16]. The last possibility can be realized by modifying the Kahler potential, or 
by adding a term ~ to the superpotential. 
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Another interesting possibility is to take the inflaton mass just a little bit smaller than 
m ~ 6 X 10~^, to decrease the amplitude of the standard inflaton perturbations. Then one 
may compensate for this decrease by adding a small contribution of the curvaton fluctuations. 
This will result in a smaller amplitude of tensor modes and a larger spectral index Ug, which 
would improve the agreement of the predictions of the simplest chaotic inflation models with 
the WMAP data. Also, as our calculations demonstrate, for certain values of parameters 
even a small contribution of the curvaton perturbations may dramatically increase the non- 
gaussianity of the combined spectrum of perturbations of metric. 

Thus, whereas the curvaton models are more complicated than the single- field inflation- 
ary models, they make the resulting scenario much more flexible, which may be important 
for a proper interpretation [22] of the coming observational data. 

Our final comment deals with the topological features of the distribution of pertur- 
bations in the curvaton scenario. We point out that in the theory of a single-component 
real curvaton field, the regions of the universe with large non-gaussianity form domain walls 
[7], reminiscent of the exponentially thick cosmic domain walls. Meanwhile in the theory 
of a complex curvaton field, which was studied in the present paper, the regions of large 
non-gaussianity form exponentially thick cosmic strings. In more complicated theories, these 
regions may form separate islands of large local non-gaussianity, resembling global monopoles. 
Since these effects have a non-perturbative, topological origin, non-gaussianity in the cur- 
vaton scenario cannot be fully described by such tools as the familiar perturbation theory 
parameters /nl and gjsih- 
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